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Abstract 
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only depend on the action functional and, therefore, unlike previously proposed 
ones, is free from any relevant ambiguity. 
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Introduction 



First order Lagrangian mechanics can be generalized to higher order Lagrangian field 
theory. Moreover, the latter has got a very elegant geometric (and homological) for- 
mulation (see, for instance, [T]) on which there is general consensus. On the other 
hand, it seems that the generalization of Hamiltonian mechanics of Lagrangian systems 
to higher order field theory presents some more problems. Several answers have been 
proposed (see, for instance, [2l|3llll[5l|6l[71|8l[9l [10] and the references therein) to the 
question: is there any reasonable, higher order, field theoretic analogue of Hamiltonian 
mechanics? In our opinion, none of them is satisfactorily natural, especially because of 
the common emergence of ambiguities due to either the arbitrary choice of a coordinate 
system [2] or the choice of a Legendre transform [TJ [HI [10] . Namely, the latter seems not 
to be uniquely definable, except in the case of first order Lagrangian field theories when 
a satisfactory Hamiltonian formulation can be presented in terms of multisymplectic 
geometry (see, for instance, [TT] - see also [12] for a recent review, and the references 
therein) . 

Nevertheless, it is still desirable to have a Hamiltonian formulation of higher order 
Lagrangian field theories enjoying the same nice properties as Hamiltonian mechanics, 
which 1) is natural, i.e., is independent of the choice of any structure other than the 
action functional, 2) gives rise to first order equations of motion, 3) takes advantage of 
the (pre-)symplectic geometry of the phase space, 4) is a natural starting point for gauge 
reduction, 5) is a natural starting point for quantization. The relationship between the 
Euler-Lagrange equations and the Hamilton equations deserves a special mention. The 
Legendre transform maps injectively solutions of the former to solutions of the latter, 
but, generically, Hamilton equations are not equivalent to Euler-Lagrange ones [TT] . 
However, the difference between the two is a pure gauge and, therefore, it is irrelevant 
from a physical point of view. 

In this paper we achieve the goal of finding a natural (in the above mentioned 
sense), geometric, higher order, field theoretic analogue of Hamiltonian mechanics of 
Lagrangian systems in two steps: first, we find a higher order, field theoretic analogue 
of the Skinner and Rusk "mixed Lagrangian-Hamiltonian" formalism [T31 [HI [U] (see 
also [IS]), which is rather straightforward (see [TT] for a different, finite dimensional ap- 
proach, to the same problem) and, second, we show that the derived theory "projects 
to a smaller space" which is naturally interpreted as phase space. Local expressions of 
the field equations on the phase space are nothing but de Bonder equations [2] and, 
therefore, are naturally interpreted as the higher order, field theoretic, coordinate free 
analogue of Hamilton equations. A central role is played in the paper by multisym- 
plectic geometry in the form of partial differential (PD, in the following) Hamiltonian 
system theory, which has been developed in [18] . 

The paper is divided into nine sections. The first four sections contain reviews of the 
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main aspects of the geometry underlying the paper. They have been included in order 
to make the paper as self-consistent as possible. The next five sections contain most of 
the original results. 

The first section summarizes the notations and conventions adopted throughout the 
paper. It also contains references to some differential geometric facts which are often 
used in the subsequent sections. Finally, in Section [T] we briefly review the Skinner- 
Rusk formalism [H] . Section [2] is a short review of the geometric theory of partial 
differential equations (PDEs) (see, for instance, [SO])- Section |3] outlines the properties 
of the main geometric structure of jet spaces and PDEs, the Cartan distribution, and 
reviews the geometric formulation of the calculus of variations [T]. Section H] reviews the 
theory of PD-Hamiltonian systems and their PD-Hamilton equations [TH]. Moreover, 
it contains examples of morphisms of PDEs coming from such theory. These examples 
are presented here for the first time. 

In Section [5] we present the higher order, field theoretic analogue of Skinner-Rusk 
mixed Lagrangian-Hamiltonian formalism for mechanics. In Section E] we also discuss 
the relationship between the field equations in the Lagrangian-Hamiltonian formalism 
(now on, ELH equations) and the Euler-Lagrange equations. In Section |6] we discuss 
some natural transformations of the ELH equations. As a byproduct, we prove that 
they are independent of the choice of a Lagrangian density, in the class of those yielding 
the same Euler-Lagrange equations, up to isomorphisms. ELH equations are, therefore, 
as natural as possible. In Section [7] we present our proposal for a Hamiltonian, higher 
order, field theory. Since we don't use any additional structure other than the ELH 
equations and the order of a Lagrangian density, we judge our theory satisfactorily natu- 
ral. Moreover, the associated field equations (HDW equations) are first order and, more 
specifically, of the PD-Hamilton kind. In Section [8] we study the relationship between 
the HDW equations and the Euler-Lagrange equations. As a byproduct, we derive a 
new (and, in our opinion, satisfactorily natural) definition of Legendre transform for 
higher order, Lagrangian field theories. It is a non-local morphism of the Euler-Lagrange 
equations into the HDW equations. Finally, in Section [9] we apply the theory to the 
KdV equation which can be derived from a second order variational principle. 

1 Notations, Conventions and the Skinner-Rusk 
Formalism 

In this section we collect notations and conventions about some general constructions 
in differential geometry that will be used in the following. 

Let be a smooth manifold. If L C A^ is a submanifold, we denote hj ii : L ^ N the 
inclusion. We denote by C°°{N) the M-algebra of smooth, M-valued functions on A^. We 
will always understand a vector field X on A^ as a derivation X : C°°{N) — y C°°{N). 
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We denote by D{N) the C°^{N)-modnle of vector fields over A^, by A(M) = @^A''{N) 
the graded M-algebra of differential forms over N and by d : A(A^) — y A{N) the 
de Rham differential. If F : A^^i — > is a smooth map of manifolds, we denote by 
F* : A(A^) — > A.{Ni) the pull-back via F. We will understand everywhere the wedge 
product A of differential forms, i.e., for uj,uji G A(A^), we will write uui instead of 
oj Aui, . 

Let a : A — y N be an affine bundle (for instance, a vector bundle) and F : A^^i — y N 
a smooth map of manifolds. Let £/ be the affine space of smooth sections of a. For 
a & ^ and x G A^ we put, sometimes, := a{x). The affine bundle on A^^i induced by 
a via F will be denoted by F°{a) : F°{A) — y N: 



F°{A) 



F°{a) 



iVi 



■^A 



and the space of its sections by F°{£/). For any section a G there exists a unique 
section, which we denote by F°{a) G F°{£/), such that the diagram 



F°{A) 

F°(a) 



-^A 



commutes. If F : A^^i — y N is the embedding of a submanifold, we also write • |i? 
(or, simply, • IatJ for F°( • ), and refer to it as the restriction of " • " to Ni (via F), 
whatever the object " • " is (an affine bundle, its total space, its space of sections or a 
section of it). 

We will always understand the sum over repeated upper-lower (multi) indexes. Our 
notations about multiindexes are the following. We will use the capital letters /, J, K 
for multiindexes. Let n be a positive integer. A multiindex of length k is a /ctuple of 
indexes I = {ii, . . . ,ik), ii, . . . ,ik < n. We identify multiindexes differing only by the 
order of the entries. If / is a multiindex of lenght k, we put \I\ := k. Let I = {ii, . . . ,ik) 
and J = (ji, . . . be multiindexes, and i an index. We denote by IJ (resp. Ii) the 
multiindex {ii, . . . ,ik, ji, . . . Ji) (resp. {ii, . . . 

We conclude this section by briefly reviewing those aspects of the Skinner-Rusk for- 
malism for mechanics fi3\ [HI [15] that survive in our generalization to higher order fleld 
theory. 

Let Q be an m-dimensional smooth manifold and q^,. . . ,q"^ coordinates on it. Let 



L G C°°{TQ) be a Lagrangian function. Consider the induced bundle Tq := Tq{tq) 



T^iT*Q) — y TQ from the cotangent bundle : T*Q — y Q to Q, via the 
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tangent bundle tq : TQ — > Q. Let q : T"!" — > T*Q be the canonical projection (see 
Diagram ([1])) 

4 b ■ 

On there is a canonical function h G C°°(T^) defined by h{v,p) := p(t>), f G TgQ, 
p e T;g, g e g. Consider also the function El := h - t^q*{L) e C°°{T^). El is locally 
given by El := Piq'' — L, where . . . ,q\ . . . ,q\ . . . ,pi, . . . are standard coordinates on 
T^. Finally, put u := g*(wo) ^ A^(T'f), ^ A^(T*g) being the canonical symplectic 
form on T*Q, which is locally given by Uq = dpidq^. w is a presymplectic form on 
whose kernel is made of vector fields over which are vertical with respect to the 
projection q. In the following, denote by I C M a generic open interval. For a curve 
7 : I 9 t I — )■ 7(t) G T\ consider equations 

^^7°(u;)+7°(rfEi) = 0, (2) 

where 7 G 7°(D(T"I')) is the tangent field to 7. Equations ([2]) read locally 




In particular, for any solution 7 of Equations (E]) as above, tq o o 7 : 1 — v Q is 
a solution of the Euler-Lagrange equations determined by L. Notice that solutions of 
Equations ([2]) can only take values in the submanifold <zT'^ defined as 

^ := {P G : there exists S G TpT^ such that i^up + {dEL)p = 0}, 

and that ^ is nothing but the graph of the Legendre transform EL : TQ — > T*Q. 
Finally, consider ^0 := q{^) C T*Q. If ^0 C T*g is a submanifold and g : ^ — ^ ^0 
a submersion with connected fibers, then there exists a (unique) function H G C°°(^o) 
such that q*{H) = El\^- Thus, for a curve cr : I 3 t 1 — > a{t) G we can consider 
equations 

Ua°iuo) + a°{dH) = 0, (3) 

where & G (T°(D(T'^)) is the tangent field to a. For any solution 7 : I — ¥ Q of the 
Euler-Lagrange equations, EL 07:! — y is a solution of Equations ([3]). If the 
map q : — > T*Q has maximum rank (which happens iff the matrix Wd"^ L / dq^ dq^H 
has maximum rank, i.e., EL is a local diffeomorphism) , then C T*g is an open 
submanifold, if is a local function on T*Q, and Equations (|3]) read locally 

_d_ j dH 

dt " (9pi 
An = — ^ ' 



5 



which are Hamilton equations. In this case, for any solution a : I — )■ T*Q of Equations 
Tq o a : I — > Q is a solution of the Euler-Lagrange equations. 



2 Geometry of Differential Equations 

In this section we recall basic facts about the geometric theory of PDEs. For more 
details see, for instance, [20] . 

Let TT : E — y M be a fiber bundle, dimM = n, dimE = m + n. In the following 
we denote hj U C M a generic open subset. For < / < A; < oo, let tta,. : J^n — > M 
be the bundle of fc-jets of local sections of vr, and tt^^i : J^vr — y J'vr the canonical 
projection. For any local section s : U — y E of it, we denote by j^s : U — y J'^tt 
its kth jet prolongation. For x E U, put [s]^ := {jks){x). Any system of adapted to 
TT coordinates (. . . , x', . . . , m", . . .) on an open subset V of E gives rise to a system of 
jet coordinates on vr^o(^) C J'^tt which we denote by (. . . , x*, . . . , . . .) or simply 
{. . . , x\ . . . , uf, . . .) if this does not lead to confusion, \I\ < k, where we put Uq := m", 
a = 1, . . . ,m. 

Now, let k < oo, tq : To — y J^vr be a vector bundle, and {. . . , x\ . . . , uf, . . . , t>", . . .) 
adapted to tq, local coordinates on Tq. A (possibly non-linear) differential operator of 
order < k 'acting on local sections of tt, with values in Tq ' (in short 'from vr to Tq ') is 
a section $ : J^n — y Tq of Tq. 

Let vr' : E' — y M be another fiber bundle and (f : E — y E' a morphism of bundles. 
For any local section s : U — y E of it, (pos : U — y E' is a local section of vr'. Therefore, 
for all < < oo, (p induces a morphism jkf : J^tt — t- J^tt' of the bundles iTk and vr^ 
defined by {jkf)[s]x '■= [v^ ° s]^, x G f/. Diagram 




commutes for all < k < I < oo. jkf is called the kth prolongation of (f. 

The above construction generalizes to differential operators as shown, for instance, 
in [22]. If $ is a differential operator as above, we denote by its Ith prolongation. 
Moreover, put := {6 G J'^tt : ^{9) = 0}. is called the (system of) PDE(s) 
determined by $. For < / < oo put also S^^ := C J^^+'vr. If (9$ is locally defined 
by 

$«(...,x\...,<,...) = 0, a = l,...,p (4) 
...,$":= $*(f'^), . . . being local functions on J'^tt, then ^^'^ is locally defined by 

{Dj^'^){...,x\...,u'},...) = 0, a = l,...,p, \J\<1, (5) 
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where D(^j-^^^,,,j^-) := Dj^ o • • -oDj^, and Dj := d/dx^ +u'}jd/du'} is the jth total derivative, 

= l,...,m. S'^^^ is called the Ith prolongation of the PDE S'^. In the 
following we put := d/duf, a = 1, . . . ,m. 

A local section s of vr is a (local) solution of ^$ iff, by definition, imj^s C ^$ or, 
which is the same, imj^+^s C S"^^^ for some / < oo. Notice that the ooth prolongation 
of (S'l^^ C J°°7r, is an inverse limit of the sequence of maps 

M < ^^i^ 4') ^^i^ 4^+^) ^ (6) 

and consists of "formal solutions" of i.e., possibly non- converging Taylor series 
fulfilling for every I. 

J°°iT is not a finite dimensional smooth manifold. However, it is a pro-finite dimen- 
sional smooth manifold. For an introduction to the geometry of pro-finite dimensional 
smooth manifolds see [21] (see also [22j, and [231 123] fo^' different approaches). In the 
following we will only consider regular PDEs, i.e., PDEs such that 4°°^ J°°7r is a 
smooth pro-finite dimensional submanifold in J°°7r, i.e., T[oo,i{(^i°°^) C J'tt is a smooth 
submanifold and tt^+i ^ : iToo,i+i{<S'i°°^) — > 'n'oo,i{'^i°°^) is a smooth bundle for all / > 0. 

There is a dual concept to the one of a pro-finite dimensional manifold, i.e., the 
concept of a filtered smooth manifold which will be used in the following. We do not 
give here a complete definition of a filtered manifold, which would take too much space. 
Rather, we will just outline it. Basically, a filtered smooth manifold is a(n equivalence 
class of) set(s) ^ together with a sequence of embeddings of closed submanifolds 

^0 ^ > ^1 ^ ■ — > ■ ■ ■ ^ > > Gk ^ > ■■■ (7) 

and inclusions i^ '■ k> 0, such that if (together with the i^'s) is a direct limit 

of (I7|) . The tower of algebra epimorphisms 

C-(^o) < < '^'^^^ C-(^,+i) < (8) 

is associated to sequence ([7]). We define C°°((^) to be the inverse limit of the tower 
([8]). Every element in C°°(i^) is naturally a function on iff. Thus, we interpret C°°{ff) 
as the algebra of smooth functions on G . Clearly, there are canonical "restriction 
homomorphisms" i\ : C°°(i^) — > C°°(^fc), A; > 0. Differential calculus over G may 
then be introduced as differential calculus over C°°(^) [2T] respecting the sequence 
([H]). Since the main constructions (smooth maps, vector fields, differential forms, jets 
and differential operators, etc.) of such calculus and their properties do not look very 
different from the analogous ones in finite-dimensional differential geometry we will 
not insist on this. Just as an instance, we report here the definition of a differential 
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form uj on ff: it is just a sequence of differential forms Uk G A{ffk), k > 0, such that 
il-i^ki^k) = uJk-i for all k. 

Finally, notice that, allowing for the i^^'s in ([7]) to be pro-finite dimensional mani- 
folds, we obtain a more general object than both a pro-finite dimensional and a filtered 
manifold. We will generically refer to such an object as an infinite dimensional smooth 
manifold or even just a smooth manifold if this does not lead to confusion. Our main ex- 
ample of such a kind of infinite dimensional manifold will be presented in the beginning 
of Section [51 



3 The Cartan Distribution and the Lagrangian For- 
malism 

Let vr : E — > M and $ be as in the previous section. In the following we will simply 
write J'' for J^vr, k < oo, and S' for S' will be referred to simply as a PDE 

(imposed on sections of vr) if this does not lead to confusion. Notice that for $ = 0, 

Recall that J°° is canonically endowed with the Cartan distribution [20] 

which is locally spanned by total derivatives, Di, i = 1, . . . ,n. ^ is a flat connection in 
TToo which we call the Cartan connection. Moreover, it restricts to cS' in the sense that 
^0 C TqS" for any 9 & S" . Therefore, the (inflnite prolongation of) any PDE is naturally 
endowed with an involutive distribution whose n-dimensional integral submanifolds are 
of the form joo-s, with s : U — >■ E a (local) solution of In the following we will 
identify the space of n-dimensional integral submanifolds of ^ and the space of local 
solutions of Sq,. 

Let tt' : E' — > M be another bundle and <f ' C J°°7r' (the inflnite prolongation of) a 
PDE imposed on sections of vr'. A smooth map F : (S' — y S is called a morphism of 
PDEs iff it respects the Cartan distributions, i.e., {dg'F){^g') = '^pie') for any 9' G S". 
The idea of non-local variables in the theory of PDEs can be formalized geometrically 
by special morphisms of PDEs called coverings [25] (see also [26]). A covering is a 
morphism ip : — > S of PDEs which is surjective and submersive. A covering 
i\} : <f — > S clearly sends local solutions of to local solutions of S . If there exists 
a covering ^ip : S — > S of PDEs we also say that the PDE S covers the PDE S' (via 
ip ). Fiber coordinates on the total space S' of a covering ip : — > ^ are naturally 
interpreted as non-local variables on S . Also notice that given a solution s of the PDE 
(o', a covering tp : S — > S determines a whole family of solutions of S "projecting 
onto s via so that ip may be interpreted, to some extent, as a flbration over the 
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space of solutions of S . Many relevant constructions in the theory of PDEs (including 
Lax pairs, Backlund transformations, etc.) are duly formalized in geometrical terms by 
using coverings. 

The Cartan distribution and the fibered structure tToo : — M of J°° determine 
a splitting of the tangent bundle TJ°° — > J°° into the Cartan or horizontal part ^ 
and the vertical (with respect to tToo) part. Accordingly, the de Rham complex of J°°, 
(A(J°°),(i), splits in the variational bi-complex {^*K ® A,d,d^), (here and in what 
follows tensor products will be always over C°°(J°°) if not otherwise specified), where 
^'A and A* are the algebras of Cartan forms and horizontal forms respectively. dX and 
d are the vertical and the horizontal de Rham differential, respectively (see, for instance, 
PU] for details). The variational bicomplex allows a cohomological formulation of the 
calculus of variations [H |20l EH [19] . In the second part of this section we briefly review 
it. 

1 1 the following we will understand isomorphism A( J°°) ~ 'if' A ® A. The complex 

> C~(J~) Ai ^ A" A'+i ■ ■ ■ 

is called the horizontal de Rham complex. An element =Sf G A" is naturally interpreted 
as a Lagrangian density and its cohomology class [^] e := i7"(A, d) as an action 
functional on sections of vr. The associated Euler-Lagrange equations can then be 
obtained as follows. 
Consider the complex 

^'^A^^-^'^Ai ® A^^^ ^^Ai® A^^^---, (9) 

and the C°°( J°°)-submodule C ^A^ (g) A" generated by elements in '^A^ (g) A" fl 
A(J^7r). x"'" is locally spanned by elements {du'^ — ufdx^) ® d"'x, where we put d"'x := 
dx^ ■ ■ ■ dx"'. 

Theorem 1 Complex Ij^ is acyclic in the qth term, for q ^ n. Moreover, for any 
u G ^A^®A" there exists a unique element E^^ G x"!" C "^A^^A^ such that E^^—uj = d-d 
for some G 't^A^^A^'^ and the correspondence i^"(^A^ ® A, d) 3 [u] i — > E^^ G is 
a vector space isomorphism. In particular, for u = d^J^, =Sf G A" being a Lagrangian 
density locally given by ^ = Ld"'x, L a local function on C°°{J°°), E{^) := E^^ is 
locally given by E{^) = j^idu"" - ufdx') ® d"x where := {-y^^Did^L are the 
Euler-Lagrange derivatives of L. 

In view of the above theorem, E{J^) does not depend on the choice of ^ in a 
cohomology class G and it is naturally interpreted as the left hand side of the 
Euler-Lagrange (EL) equations determined by =Sf . In the following we will denote by 
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S'el C J°° the (infinite prolongation of tlie) EL equations determined by a Lagrangian 
density. Any ■& e "^A^ ® A""^ sucli tliat 

E{^) - rf^if = M (10) 

will be called a Legendre form (TU]. Equation f lTU]) may be interpreted as the first 
variation formula for the Lagrangian density =Sf. In this respect, the existence of a 
global Legendre form was first discussed in [ 27] . 

Remark 1 Notice that, if d & '^h}®K^~^ is a Legendre form for a Lagrangian density 
G A", then °d + c?^^ is a Legendre form for the d-cohomologous Lagrangian density 
S£ + dq, Q G A"~^, which determines the same EL equations as S£ . Moreover, any 
two Legendre forms i), d' for the same Lagrangian density differ by a d-closed, and, 
therefore, d-exact form, i.e., d — d' = dX, for some A G ® A^~'^. 

Remark 2 Finally, notice that complex ^ restricts to holonomic sections jooS ofn^, 
s being a local sections of tt, in the sense that, for any such s, there is a (unique) 
complex 

> ^A^l, ^ ^A^ ® Ai I, > ^A^ ® A^j (11) 

where j := joo-s, such that the restriction map 'tfA^ 0A — > ^A^ ^Mj — ^-^^Ij ®c°°{m) 
A(M) is a morphism of complexes. Moreover, complex / flT]) is acyclic in the qth term 
and the correspondence defined by H"-(jifA^ ^ A"'\j,d\j) 3 [uj\j] i — > G oj G 

^A^ ® A", is a vector space isomorphism. 

4 Partial Differential Hamiltonian Systems 

In [18] we defined a PD analogue of the concept of Hamiltonian system on an abstract 
symplectic manifold which we called a PD- Hamiltonian system. In this section we 
briefly review those definitions and results in p8] which we will need in the following. 

Let a : P — > M be a fiber bundle, A := C°°{P), x^,...,x^ coordinates on M, 
dimM = n, and g^, . . . , g*" fiber coordinates on P, dimP = n + m. Denote by C{P, a) 
the space of (Ehresmann) connections in a. C{P, a) identifies canonically with the space 
of sections of the first jet bundle ai^o '■ J^ot — > P and in the following we will understand 
such identification. In particular, for V G C(P, a), we put . . . , := V*(g^),..., 
. . . , qf^ . . . being jet coordinates on J^a. 

Denote by Ai = ^^A'l C A(P) the differential (graded) ideal in A(P) made of 
differential forms on P vanishing when pulled-back to fibers of a, by Ap = @f^A^ its 
p-th exterior power, p > 0, and by VA(P, a) = 0^VA'^(P, a) the quotient differential 
algebra A(P)/Ai, : VA{P,a) — )■ VA{P,a) being its (quotient) differential. 
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Remark 3 For instance, if a = vToo : P = — > M , then, using the Cartan 
connection ^ G C{J°° ,tToo) , one can canonically identify VA^{J°^ ^Hoo) with "loA^ 
and with the vertical de Rham differential. More generally, for any k > 0, 
yA^( J^, vrjfc) 0c'^(^jk^-j C°°(J^+V) identifies canonically with the C°°{J'^^^'n) -module 
n A(J^+V) of [k + l)th order Cartan forms. 

Now, for any A; > 0, put fi^(P,a) := h^tT^ and n^{P,a) := fi'=(P, a)/A^+"-i. It is 
easy to show that V^{P, a) ~ VK^{P, a) ®a Ki-i- An element oo G i7^(P, a) determines 
an affine map 

C(P, a)3V ^ i^u := pv(u;) G VA'=-1(P, a) ®a (12) 

where 

Pv:A(P)^yA^-i(P,a) ®a K 

is the canonical projection determined by the connection V. The linear part of the 
affine map ([T^ naturally identify with the class u + A^+"-^ in ^t{P, a) (see [IB] for 
details). Notice that, since f lT2|) is affine, it is actually point- wise and, therefore, can be 
restricted to maps. Namely, if P : Pi — > P is a smooth map, □ G P°(C(P, a)), then 
an element iuF°{uj) G P°(VA^^^(P, a) ®a AJ^) is defined in an obvious way. 

Definition 1 A PD-Hamiltonian system on the fiber bundle a : P — > M is an element 
uj G fi^(P, a) such that duo = The first order PDEs 

^iicr^lcr = 

on (local) sections a of a are called the PD-Hamilton equations determined by u. 
Geometrically, they correspond to the submanifold 

.= g ji^t : i^ujp = 0,p = aifl{9)} C J^a. 

Let w be a PD-Hamiltonian system on the bundle a : P — > M and consider the subset 
Pi := ai o(<Cj°^) C p. In the following we will assume Pi C P to be a submanifold and 
Oi := q;|pj : Pi — > M to be a subbundle of a. ai is called the first constraint subbundle 

of UJ. 

As an example, consider the following canonical constructions. Let a : P — > M 
be a fiber bundle and ...,q^,... fiber coordinates on P. VL^{P,a) (resp. il^(P, «)) 
is the C°°(P)-module of sections of a vector bundle ^^a : — > P (resp. r^a : 
J'^a — > P), called the multimomentum bundle of a (resp. the reduced multimomentum 
bundle of a). Recall that there is a tautological element 9^ G Q}{.y£a, iia) (resp. 
Gq, G V^^J^a.r'^a)), where /xa := a o /igo; (resp. T^a := a o Tga), which in standard 
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coordinates . . . ,x\ . . . , q^, . . . , p\, . . . ,p on (resp. . . . ,x\ . . . , q^, . . . , p\, ... on 

J'^a) is given by 

0Q = p\dq^(r~'^Xi - pcTx (resp. 6q = p\d^q^ ® cT'^Xi), 

where d"'~^Xi := iQ/Qxid"'x [28\. dOa is a PD-Hamiltonian system on fia locally given 
by 

dOa = dp^dq'^d'^ ^Xi - dpd^x. 
Notice that dOa determines empty PD-Hamilton equations. 

Example 1 A PD-Hamiltonian system is canonically determined, on the fiber bundle 
a : P — y M , by the following data: a connection V in a and a differential form 
=^ G AJ^. Let . . . , g"^, ... be fiber coordinates in P and . . . ,x\ . . . , q^, . . . ,p\, . . . ,p 
(resp. . . . ,x\ . . . ,q^, . . . ,p\, . . .) standard coordinates in .Ma (resp. J^a). Let ^ 
he locally given by ^ = Ld^x, L a local function on P. Obviously, V determines a 
section Ey : J'^o. — of the projection — > J'^a, which in local standard 
coordinates reads = p\Vf. Put 6v := Sy(0Q,). In local standard coordinates, 

9v = p\dq^d"~^Xi — p\Vfd'^x. Put also. 

Locally, 6^^v = p^dq^d^'^Xi — Ej^^yd^x, where E^y := p\Vf — L. Finally, consider 
|^jsf,v '■= dO^y. Locally, 

|^if,v = dp^dq'^d'^'^Xi - dE^^^d'^x. 

^y.v is the PD-Hamiltonian system on r^a determined by V and The associated 
PD-Hamilton equations read locally 

/ PAn = - PB^^f 

where we denoted by "•ji " the partial derivative of " with respect to the ith independent 
variable x^ , i = 1, . . . ,n. 

We conclude this section by discussing two examples of morphisms of PDEs coming 
from the theory of PD-Hamiltonian systems. 

Example 2 Let a : P — > M be a fiber bundle, u a PD-Hamiltonian system on it, 
a' : P' — M another fiber bundle, (3 : P' — > P a surjective, submersive, fiber bundle 
morphism, andu' := P*{uj). u' is a PD-Hamiltonian system on a' . Denote by S C J°°a 
(resp. S" C J^a') the ooth prolongation of the PD-Hamilton equations determined by 
(jj (resp. uj' ). We want to compare S and S" . In order to do this, notice, preliminarily, 
that J°°a' covers J°°a via joo/3 : J°°a' — )■ J°°a. Moreover, it can be easily checked 
that a local section a' of a' is a solution of S" iff the section [3 o a' of a is a solution of 
S . We now prove the formal version of this fact. 
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Proposition 2 {jooP){^') C S and jooP '■ (S" — > £' is a covering. 

Proof. Consider j 1 13 : J^a' — > J^a. It is easy to check that (fj?^ = C 
J^a'. Similarly, S" = (joo/3)^^(<^) C J°°a'. In particular, j^P : S" — )■ ^ is the 
^^restriction" of j^oP '■ J°°a' — > J°°a to £' C. J'^a and, therefore, is a covering. ■ 

Example 3 Let a : P — > M , u and S C J°°« he as in the above example, and 
a\ : Pi — > M the first constraint suhhundle of u. Put Ui := i*p_^{u). Ui is a PD- 
Hamiltonian system on a\. Denote by S\ C J°°ai the ooth prolongation of the PD- 
Hamilton equations determined by ui . We want to compare S and Si . In order to do 
this, notice, preliminarily, that J'^ai may be understood as a submanifold in J°°a via 
jooipi '■ J°°ai ^ J°°a. Moreover, it can be easily checked that any solution of S is also 
a solution of Si (while the vice-versa is generically untrue). We now prove the formal 
version of this fact. 

Proposition 3 S G Si. 

Proof. Recall that the projection ai^Q : J^a — )■ P sends SS^^ to Pi. As a consequence, 
S C J°°ai. Moreover, by definition of ooth prolongation of a PDE, it is easy to check 
that 

S = Sf} J°°ai 

= {9 = [a]- G J^ai 

= Si. 



5 Lagrangian-Hamiltonian Formalism 

In this section we show that the Skinner- Rusk mixed Lagrangian-Hamiltonian formalism 
for first order mechanics [IHl HU [15] (see Section [1]) is straightforwardly generalized to 
higher order Lagrangian field theories. 

First of all, let us present our main example of a filtered manifold. Let vr : E — > M 
be a fiber bundle. Consider the infinite jet bundle ir^ : J°° — )■ M for which = A"^, 
g > 0. Moreover, the C°°( J°°)-module fi^(J°°,7roo) ^ ^A^^A""! is canonically filtered 
by vector subspaces Wk := '^A^ ® A"-^ n A( J^^+^vr), /c > 0. Denote by fi^ C J°°, tToo) 
the J°°)-submodule generated by Wk-, k > 0. Then, for all k, is canonically 
isomorphic to C°°{J°°) ^fc 

^llcnlc--- cnlc nl^i c ■ • • c n\J°°, ttoo), (i3) 



[inM<r]T = 0, xeM} 
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is a sequence of C°°( J°°)-submodules. Notice that, for any k, Ofc is the module of 
sections of a finite-dimensional vector bundle — > J°°. Moreover, the inclusions (fT3il 
determine inclusions 

Jo^ C 4 C ■ ■ ■ C 4 C C ■ • • 

of vector bundles. := IJ^ 4 is then an infinite dimensional (filtered) manifold and 
the canonical projection Tq : J"^ — > J°° an infinite dimensional vector bundle over 
whose module of sections identifies naturally with fi^( J°°, tToo). We conclude that 
Tq : — > J°° is naturally interpreted as the reduced multimomentum bundle of tToo- 
Denote by . . . , x*, . . . , u", . . . . . . standard coordinates on J'^ . We will also consider 

the bundle structures 4 — ^ ^) ^ > 0, and r"!" := tToo o : J'^ — )■ M. 

In the following we denote by U' C J°° a generic open subset. Notice that any 
(local) element G 'ioA^ A"~^ = f2^(J°°, vToo), in particular a (local) Legendre form, is 
naturally interpreted as a section § : U' — )■ J"'' of Tq. Put then . . . , -i?^* := ^*{Pa^), ■ ■ ■ 
which are local functions on J°° such that 'd = 'd^^'^{duf — uj^dx^) ® It follows 

that, locally, 

M = -{Di^l;' + S^ji^^J){du'} - u%dx') ® d'^x e ^A^ ® A", 

where 5^ = if / ^ K, while 5^ = 1 if / = if . 

Now, in Example [H put a = tToo '■ P = J°° — > M and V = the Cartan connection 
in TToo. ^ G A" = A" is then a Lagrangian density in tt. Put S := S<^, 6^ := Q^^^ and 
:= cjj^','^. cjjf is a PD-Hamiltonian system on r'^ : — )■ M canonically determined 
by ^ . Locally, 

UJ^ = dpl;''du'}d'^~^Xi - dEyd^x, 

where := p^^uj.- — L. Let a : U — .P be a local section of t\ and j := o a : 
U — > J^. Put ... , af := a*{uf) = j*{uf), . . . , a^J := a*(p^-*), . . . which are local 
functions on M. Then, locally, 

and the PD-Hamilton equations determined by read locally 

r p'^\ = diL-5'j,pi^ _ 

Uj ,i= 

We call such equations the Euler-Lagrange-Hamilton (ELH) equations determined by 
the Lagrangian density ^ . Notice that they are first order PDEs (with an infinite 
number of dependent variables). Denote by Selh C J°°r^ their infinite prolongation. 
In the following theorem we characterize solutions of S'elh- As a byproduct, we derive 
the relationship between the ELH equations and the EL equations. 
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Theorem 4 A local section a : U — f .P of is a solution of the ELH equations 
determined by the Lagrangian density S£ iff it is locally of the form a = -d o j°°s where 
1) s : U — > E is a solution of the EL equations S'el and 2) : U' — )■ P is a Legendre 
form for =Sf . 

Proof. Let a : U — > be a local section of First of all, let cr be of the form 
a = d o j where 1) j : U — > J°° is a local section of tToo and 2) d -.U' — )■ is a local 
section of Tq : P — > J°° . Then, 

Therefore, locally, 

where ■ ■ ■ ,jf '■= j*{uf), . . . and they are local functions on M. Thus, if ■(9 is a Legendre 
form and j = jooS for some local solution s : U — > E of the EL equations then, in 
particular, = jf^, a = 1, . . . ,m, i = 1, . . . ,n, and 

On the other hand, let a : U — > P be a local section of r''" and j := Tq o a : f/ — y J°°. 
Locally, there always exists a section d : U' — )■ P of Tq, such that a = 'd o j. Notice, 
preliminarily, that d is not uniquely determined by a except for its restriction to imj. 
If 0" is a solution of the ELH equations then, locally, 

= = (d^ + d''^)\, + (jT,. -3%)dVc.\ ® d'^x. 

Since {d^Pa^)\a ® d"'x and {d{} + d^^)\j are linearly independent, it follows that 

r (dd + d^^)\j = 

t Jin = ill- 

In particular, j = jooS, where s = 71^0,0 ° j- 

Now, let "i^o be a Legendre form for =Sf . Then d^^ = E{^) — d^Q and, therefore, 
{d'd — d'do + E{^))\j = 0. Recall that d restricts to j = jooS (Remark |2]). Thus, 

d\,{^-^o)\j = E{^)\,. 

In particular, E[^) \j is (i|j-exact. In view of Remark[2l this is only possible if E{^) \j = 
0, i.e., s is a solution of the EL equations. We conclude that 
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i.e., — '&o)\j is d\j-closed. Again in view of Remark [21 this shows that, locally, 

(^-^o)|j = d\ju\j = Mj 

for some u G 'ifA^ ® A""^. In particular, we can put = "^o + du and, therefore, ^9 is a 
Legendre form for ^ as well. ■ 
We now prove a formal version of the above theorem. Put P '■= o ° Tq : J°°t'' — > 

Theorem 5 p{<Selh) C S'el and p : S'elh — > ^el is a covering of PDEs. 

Proof. In J°°r"'' consider the submanifold (S'l made of ooth jets of (local) sections 
a : U — )■ J"!" of the form a = ^ o j^s, where s : U — )■ is a local section of vr, 
and : U' — > is a local Legendre form. It can be easily checked that S'l is locally 
defined by 

Pi;'\K^ + ^^j^Ak = DKidiL) - 6^oDk^ .... 

Clearly, the Cartan distribution restricts to S'l and, therefore, S'l can be interpreted as 
a PDE. Moreover, it is easily seen from (IT^ that S'l covers J°° via p. Denote by 

P)' d I „ a d I I.i d 

the jth total derivative on J°°t'^ , j = 1, . . . ,n. (S'elh is locally defined by 

pi;^\K^ = D'^{diL)-6^j,p'^^\K 



which is equivalent to 



p'J\K^ = DK{^iL)-5^Jy^\K 



Moreover, on S'elh 

and, therefore, Dkj^ = 0, a = 1, . . . , m. It then follows from f[T^ . that Selh = 
Sl ^P'^i^EL)- In particular, p : S'elh — >■ Sel is the "restriction" of p : c^i^ — > J°° to 
Sel C. J°° and, therefore, is a covering. ■ 
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6 Natural Transformations of the Euler- Lagrange- 
Hamilton Equations 

Properties of Legendre forms discussed in Remark [T] correspond to specific properties 
of the ELH equations which we discuss in this section. 

First of all, notice that the ELH equations are canonically associated to a Lagrangian 
density. But, how do the ELH equations change when changing the Lagrangian density 
into a (i-cohomology class? In particular, does an action functional uniquely determine 
a system of ELH equations or not? In order to answer these questions consider e 
® A"~^. {} determines an automorphism : — of the fiber bundle Tq via 

In particular, o \['^ = rj. Clearly, ^f^^ = \E'_,9. 
Lemma 6 ^^(c.;^) = - r^*{di^). 
Proof. Compute, 

= d[(^;oS*)(e) + (ro^ovI/,)*(^)] 
= rf[(vi/*oS*)(e) + rt*(^)]. 

Now, since, locally, . . . , "^^{Pa^) = Pa^ — '&a\ . . ., we have 

{%o^*){p) = {p'a-€>n- 

Thus, locally 

o E*)(e) = {p'J - - {p'J - d'^')u%d^x 

= S*(e)-ro^*(^). 

We conclude that 

%{uj,) = d[{%o^*){Q) + rl*{^)] 

= rf[S*(0)-rt*O?)+rt*(^)] 

= Uj^-Tl*{d^). 
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Theorem 7 Let =Sf' = =Sf + dg, G A" ^ , be another Lagrangian density (thus, =Sf' 
determines the same EL equations as S^). Then '^*^vj^uj^) =00^1. 

Proof. Notice, preliminarily, that 

Tl*{dd^Q) =Tl*{dd^Q) 

= —TQ*{ddg) 
= —dTQ*{dg). 

Therefore, in view of the above lemma, 

Kv.i^y) = - 4*idd^Q) 

= d[J:*{Q) + T^o*i^)]+dTl,*(dg) 

= d[i:*ie) + 4*i^ + dg)] 

= dQ _sf/ 



Corollary 8 An action [=Sf] G H", =Sf G A", uniquely determines a system of ELH 
equations, modulo isomorphisms of PD-Hamiltonian systems. 

We conclude that the ELH equations are basically determined by the sole action 
functional and not a specific Lagrangian density. 

Theorem 9 Let d G "ioh} ® K"-^^ he d-closed, hence d-exact. Then, for every La- 
grangian density ^ G A", \E'.^ is a symmetry of the ELH equations determined by ^ in 
the sense that jao^^ : J°^r'^ — )■ J°°t'^ preserves (Selh- 

Proof. By definition of infinite prolongations of a PDE and infinite prolongation of 
a morphism of bundles, it is enough to prove that ji'^^ : J^r^ — )■ J^t^ preserves 
^ELH '■~ ^ J^t'^ . Notice, preliminarily, that, in view of the proof of Theorem IH 
we have 

Now, let c G (S'^Ihj P '■= '^loi'^) ^-^^ ^ G Tpj'^ be a tangent vector, vertical with respect 
to rt. Consider also c' := (ji*^)(c), P' := ^^(P) = r|o(c') and ^ := In 
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particular, ^' G Tpij"^ is vertical with respect to r"!" as well. Let us prove that c' G (4lh- 
In view of Lemma El 

miiujj^) =u^- tI\M) =uj^- rl*{(fd). 

Compute 

VV(w^)p' = i£_ic'^l{uj^)p = i^ic{uj^)p - i^ic[r^*id^'&)]p = -i^,>i.g^{d^^)e = 0, 

where 6 = t^{P) G J°° and ^" = {dT^){0 e TgJ^^ is a tangent vector, vertical with 
respect to vToo- It follows from the arbitrariness of that ic'iojj^)p' = 0. ■ 



7 Hamiltonian Formalism 

In this section we present our proposal of an Hamiltonian formalism for higher order 
Lagrangian field theories. Such proposal is free from ambiguities in that it depends 
only on the choice of a Lagrangian density and its order. Moreover, d-cohomologous 
Lagrangians of the same order determine equivalent "Hamiltonian theories" . 

First of all, we define a "finite dimensional version" of the ELH equations (see also 
|17]). In order to do this, notice that, in view of Remark [3l for all k > 0, Wk 
is canonically isomorphic to the C°°( J'^"'"^)-module of sections of the induced bundle 
"^fc+i fc('^^^fc) — > J^~^^. We conclude that — )■ J°° is canonically isomorphic to the 
pull-back bundle tt^ ki^'^'^k) — ^ ^ k > 0. Notice that the coordinates . . . ,Pa\ ■ ■ ■, 
\I\ < k, on jI identify with the pull-backs of the corresponding natural coordinates on 
J'^TTk which we again denote by . . . . . .. 

Now, let =Sf G A" be a Lagrangian density of order / -|- 1, i.e., =Sf G A" fl A(J'+^). Let 
u'l be the pull-back of uj^ onto jj . u[ is a PD-Hamiltonian system on j/ — t- M, and 
it is locally given by 

uj[ = dpl;'du'}d''-^Xi - dEid'^x, 
\i\<i 

where Ei = X]|/|</Pa''"/j — is the restriction of to j/. Notice that ul is also 
the pull-back via jj — > T^i^niJ'^T^i) of a (unique) PD-Hamiltonian system ui on 
'^i+ii{J'^'^i) — ^ is locally given by the same formula as u'l and it is a con- 

strained PD-Hamiltonian system, i.e., its first constraint bundle — > M is a proper 
subbundle of 7ii_^-^ i^J'^ tti) — y M. Let us compute it. Let P G 7r^°+i ^(^^vr;) and 

■= '^■l_^^ll{T^nl){P) G J'+^ Then P G ^ iff there exists c in the first jet bun- 
dle of T^i+i^iiJ^T^i) — > M such that ic{^^i)p = 0, i.e., iff there exist real numbers 
. . . , Cj.j, . . . , c^'*.j, . . ., |/| < /, such that 

r ci;\. = {diLm-6^j,P;^-\ \I\ < I + 1 
lc}.. = Py„ \J\<1 
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where we put = for |/| = / + 1, and . . . , P% := u%{P), P^-' := p^-\P), . . ., 
\J\, \K\ < I, a = 1, . . . ,m. Thus, for \I\ = I + 1, P should be a solution of the system 

diL - 5'j,pt' = 0, \I\ = l + l. (16) 

Equations (|T6|) define ^ locally. 

Remark 4 is a submanifold of tTi^hIJ'^tii) of the same dimension as J^vr^, and 
— )• J^^^ an affine subbundle of 71°^^ i{J'^Tri) — > J^^^ . 

Let be the image of 0^ under the projection T^iJ^niJ^T^i) — >■ J^t^i- 

Assumption 1 We assume 0^^ to he a submanifold of J'^tti and t^t^i\^q : l^o — )■ M 
to he a smooth subbundle of t'^tti. We also assume that the projection q : ^ — > is 
a smooth submersion with connected fibers. 

Notice that, as usual, all the above regularity conditions are true if we restrict all the 
involved maps to suitable open subsets. 
The following commutative diagram summarizes the above described picture: 

< ' J} 





M 



Theorem 10 Under the regularity Assumption Ul there exists a unique PD- 
Hamiltonian system uq on — > M , such that i*^{ui) is the pull-back of via 
q: ^ — ^ ^0- 
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Proof. Since q : — > has connected fibers and i*^{oo) is a closed form, it is 
enougli to prove tiiat iyi^(a;i) = for all vector fields Y e D{j^) vertical with respect 
to q. Let Y e D(7r[_^_^^;( J'^tt;)) be vertical with respect to 7rf^_i^;( J'^'tt/) — > J^^i, and 
y := Then Y is locally of the form 

\K\=l+l 

for some . . . , y|, . . . local functions on ^. Now Y e D(^) iff, locally, 

\i\=i+i 

Compute 

\K\=l+l 111=1+1 

Thus Ei\^ is the pull-back via g of a (unique) local function H on ^q. Moreover, 

= -Y{Ei\^)(r'x = 0. 

It follows from the arbitrariness of Y that is the pull-back via q of the PD- 

Hamiltonian system loq on — > M locally defined as 

\n<i 



Definition 2 ujq is called the PD-Hamiltonian system determined by the {I + l)th 

order Lagrangian density =Sf , and the corresponding PD-Hamilton equations are the 
Hamilton-de Bonder- Weyl (HDW) equations determined by ^ . 

Definition 3 A Lagrangian density of order I is regular at the order I + 1 iff 
the map ^ — >■ J'^tt; has maximum rank. 

The Lagrangian density ^ of order / + 1 is regular at the order / -|- 1 iff the matrix 

H(L)(^) \\{d^diL){e)\\ |7|, \K\^l + 1, 

where the pairs (a, /) and K) are understood as single indexes, has maximum rank 
at every point 9 e J'+^. In its turn, this impfies that is an open submanifold of J'^tt/ 
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and, in view of Remark H] and Assumption ^ q : ^ 



^0 is a diffeomorphism. In 



particular, ujq is a PD-Hamiltonian system on an open subbundle of r'vr; locally given 



\i\<i 

where, now, H is a. local function on J^tti. In this expected, the HDW equations 

read locally 



Notice that the HDW equations are canonically associated to a Lagrangian density and 
its order and no additional structure is required to define them. Moreover, in view of 
Theorem [71 two Lagrangian densities of the same order determining the same system 
of EL equations, also determine equivalent HDW equations. Finally, to write down the 
HDW equations there is no need of a distinguished Legendre transform. Actually, the 
emergence of ambiguities in all Hamiltonian formalisms for higher order field theories 
proposed in the literature seems to rely on the common attempt to define first a higher 
order analogue of the Legendre transform and, only thereafter, the "Hamiltonian the- 
ory" . In the next section we present our own point of view on the Legendre transform 
in higher order Lagrangian field theories. 

8 The Legendre Transform 

Keeping the same notations as in the previous section, denote by ^S'elh the infinite 
prolongation of the PD-Hamilton equations determined by uji and by p' : 7?^°^^ ^ 
E the natural projection. 

Proposition 11 {jooP')i!^ELH) C S'el and 

jooV' '■ ^S'elh — ^ ^EL is a covering. 

Proof. The proof is the finite dimensional version of the proof of Theorem [5] and will 
be omitted (see also [E]). ■ 

Denote also by S"^ the infinite prolongation of the PD-Hamilton equations determined 
by and by i^j^ the infinite prolongation of the HDW equations. 

Proposition 12 (joo?) (<^if) C i^j^ and jooQ '■ — ^ ^ covering. 

Proof. It immediately follows from Theorem [10] and Proposition [2j ■ 



by 
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Notice that, in view of Propositions [3l [TT] and [121 there is a diagram of morphisms of 
PDEs, 



e>ELH ^ 



JooQ 



(17) 



whose vertical arrows are coverings. Therefore, the inclusion ^S'elh C S"^ may be 
understood as a non local morphism of (Sel into (Sh- We interpret such morphism as 
Legendre transform according to the following 



Definition 4 We call diagram [T)y the Legendre transform determined by the La- 
grangian density =Sf . 

Any Legendre form of order /, {} : J°° — > jj — > 7ii_^^ i^j'^ni), determines a section 
Joo^Uel '■ ^EL — > ^S'elh of the covering joop' : ^(Selh — > ^el and, therefore, via 
composition with joo?, a concrete map Sel — > <Sh- Nevertheless, among these maps, 
there is no distinguished one. 

We now prove that, if ^ is regular at the order / + 1, then S'h itself covers Sel- This 
result should be interpreted as the higher order analogue of the theorem stating the 
equivalence of EL equations and HDW equations for first order theories with regular 
Lagrangian (see, for instance, Let us first prove the following 

Lemma 13 If is regular at the order I + 1, then ^S'elh = <S'h ■ 

Proof. The proof is in local coordinates. Let a : U — )■ 7r^+i;(jWz) be a local section 
of T^i+iiiJ"^ T^i) — > M. Suppose imcr C Then, locally, 

diLoa- 6]^ai^ = 0, |/| = / + 1. 
Now, ij^„uji\a is locally given by 



As already outlined, the annihilator of D(i^) in A^(7r^°^]^ ^(J'I'tt;))! <^ is locally spanned 
by 1-forms 

\i:=d{diL-6yt^)\^, \I\=l + l. 
Therefore, ij^o-z^(ci;z)|o- = iff, locally, 

^hM^= Yl //AaU®^^"a;, (18) 

\i\=i+i 
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for some local functions . . . , , ... on im a, where 

:= d\diL - 5yt%^ = Y.\K\<iAdfdiLd%'^ - syVa^y, \I\ = l + 1. 

Equations (|T8|) read 

where /[i] is the number of times the index i appears in the multiindex /. Since the 
vertical forms . . . , £^uj\a, ■ ■ ■ , (i^p^'*|o-, . . . are linearly independent, ijj^ai*^{ui)\o- = iff, 
locally, 

-a^\, +diL o a - - Y.\k\=i^^ f^Kd^iL o a = 0, |/| < / + 1 

^?,-^/"=0, |/|</ , (19) 

^?n-cT% + ^m = 0, \I\=l 

for some . . . , /°, . . .. It follows from the third of Equations f[T^ that 

//. = -7FfTTK«-^/")' 1^1 = ^- (20) 
Moreover, since im cr C the first equation, for |/| = / + 1, gives 

\K\=l+l \J\=l \J\=l 

and, in view of the regularity of =Sf and Equations ( 12U]) . 

= //" = 0, |/|=/. 

Substituting again into fll9p . we finally find that the PD-Hamilton equations 
= are locally equivalent to equations 



!^\,= d[L-6Wj. |/|</ + 1 

1/1 </ 



Pa n = O^L - Ojjp- 
Uf,i = M/j, 



which are the PD-Hamilton equations determined by w/. ■ 

Now, suppose that ^ is regular at the order / + 1. Then, as already mentioned in the 
previous section, q : ^ — )■ is a diffeomorphism, and q*{u!Q) = i*^{u!i). Therefore, 
jooq : (ojf — )■ (OH is an isomorphism of PDEs and the Legendre transform f[T71) reduces 



to 



'S'ELH 



Sel <S'h 
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Moreover, maps to E via tti^ o TqTT/ and such map is a morphism of bundles 

(over M). The induced morphism J°°TqTTi — > J°° restricts to a morphism of PDEs, 
K : S'h — > locally defined as k*{u'^) = Uq\k, \K\ > 0. It is easy to show that 
diagram 

&ELH '^H 

joc.q , 

Sel ' > 

commutes, so that k = jooP' ° UooQI^h)'^ ■ Consequently, ^{S'h) C S'el and k : S'u — >• 
S'el is, a covering. Summarizing, we have have proved the following 

Theorem 14 If ^ is regular at the order I + 1, then S'h covers S'el- 

Finally, it should be mentioned that in most cases, even if the Lagrangian density is 
not regular, S^ covers Sel via n and, therefore, S"^ itself covers Sel (see the example 
in the next section). 



9 An Example: The Korteweg-de Vries Action 

The celebrated Korteweg-de Vries (KdV) equation 

0t - 600a; + 0^^^ = (21) 

can be derived from a variational principle as follows. Introduce the "potential" u by 
putting Ux = 0. Equation ([5T]) becomes the fourth order non-linear equation 

for sections of the trivial bundle vr : x M 9 (t, x; u) i — y (t, x) G M^. In its turn, fl22l) 
is the EL equation determined by the action functional 

j (ul - \uxUt + \ul,^dtdx. 

Choose the second order Lagrangian density 

^ = ("x - k^xUt + \ul^)dtdx. (23) 

Since the matrix 

/ 
H(L) =000 
\ 1 
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has rank 1, =Sf is not regular. Let t, x, u, Ut, Ux, p'*, p'^, ■ ■ ■ p^'\ ... be natural coordinates 
on J^i, i,j = t,x. Then 

Ui = dp'^dudx — dp'^dudt + dp^'^dufdx — dp^'^dutdt + dp^'^duxdx — dp^'^duxdt — dE^dydtdx 
where 

^KdV := P'\ + P'^'Ux + p^-^Utt + (p*-^ + + p"'-"'Uxx -ul + \UxUt - \u\^. 



Accordingly, ^S'elh reads 



' p-\t+p-'',x = 



2 "a; 



= 3ul - 



i,j = t,x 



(24) 



p*-* = 

pt.X I pX .t Q 

^ p''-'' -Uxx = 

which clearly cover fl22|) . Notice that the last three equations in (12^ define Thus, 
is coordinatized by t, x, u, Ut, Ux, Uu, Utx, P*'^, P^'^ and 

i*gg{uji) = dp'^dudx — dp^dudt — dp^'^{dutdt + duxdx) — dp^'^duxdt — dEi<i(iY\gsdtdx, 

where 



^Kdvl,^ := P'^ut + p-'^Ux + \{p''-''f - ul + \UxUt. 



Accordingly, reads 



f p.t^^ = 

f r 1 
P*'^,t 



-3ul + \ut 



i = t,X 



, Ux,x = P 

Notice that, even if the Lagrangian density is not regular, and variables Utt,Utx are 
undetermined, covers (!22|) . Finally, ^ is defined by the sixth and the seventh 
equations in and, therefore, it is coordinatized by t,x,u,Ut,Ux,p^'^,p^'^- Thus, Uq 
and (^HDW are given by exactly the same coordinate formulas as i*^{coi) and S'^ . In 
particular, S'hdw itself covers coel- 

Finally, recall that the KdV equation is Hamiltonian, i.e., it can be presented in the 
form Ut = A[E[^)), where is a top horizontal form in the infinite jet space of the 



26 



bundle 3 {x; u) i — )• x G M, and A is a Hamiltonian ^-differential operator (see, 
for instance, [20] )• Since Hamiltonian PDEs play a prominent role in the theory of 
integrable systems, it is worth to mention that such property (which is based on a 1+1, 
"covariance breaking" splitting of the space of independent variables (t, x)) is directly 
related with the present covariant Hamiltonian formalism as shown, for instance, in 
PU] . There the author provides a multisymplectic framework for the KdV equation by 
choosing, along the lines of [6], a "quasi-symmetric" Cartan form for the Lagrangian 
density fl23l) . Such Cartan form is unique for a second order theory. Therefore, the 
formalism of [29] is actually equivalent to ours, in the special case of a second order 
theory. 

Conclusions 

In this paper, using the geometric theory of PDEs, we solved the long standing problem 
of finding a reasonably natural, higher order, field theoretic analogue of Hamiltonian 
mechanics of Lagrangian systems. By naturality we mean dependence on no struc- 
ture other than the action functional. We achieved our goal in two steps. First we 
found a higher order, field theoretic analogue of the Skinner-Rusk mixed Lagrangian- 
Hamiltonian formalism [131 O IIS] and, second, we showed that such theory projects 
naturally to a PD-Hamiltonian system on a smaller space. The obtained Hamiltonian 
field equations enjoy the following nice properties: 1) they are first order, 2) there is a 
canonical, non-local embedding of the Euler-Lagrange equations into them, and 3) for 
regular Lagrangian theories, they cover the Euler-Lagrange equations. Moreover, for 
regular Lagrangian theories, the coordinate expressions of the obtained field equations 
are nothing but the de Bonder higher order field equations. This proves that our theory 
is truly the coordinate-free formulation of de Bonder one [2]. 
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